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ABSTRACT 

We know now from radial-velocity surveys and transit space missions that planets only a few times more massive than our Earth 
are frequent around solar-type stars. Fundamental questions about their formation history, physical properties, internal structure, and 
atmosphere composition are, however, still to be solved. We present here the detection of a system of four low-mass planets around the 
bright (V=5.5) and close-by (6.5 pc) star HD 219134. This is the first result of the Rocky Planet Search program with HARPS-N on 
the Telescopio Nazionale Galileo in La Palma. The inner planet orbits the star in 3.0937 ± 0.0004 days, on a quasi-circular orbit with 
a semi-major axis of 0.0382 ±0.0003 AU. Spilzer observations allowed us to detect the transit of the planet in front of the star making 
HD219134b the nearest known transiting planet to date. From the amplitude of the radial-velocity variation (2.33 ± 0.24ms“*) and 
observed depth of the transit (359 ± 38 ppm), the planet mass and radius are estimated to be 4.46 ± 0.47 Mg and 1.606 ± 0.086 R®, 
leading to a mean density of 5.89 ± L17gcm“^, suggesting a rocky composition. One additional planet with minimum mass of 
2.67 ± 0.59 Mg moves on a close-in, quasi-circular orbit with a period of 6.765±0.005 days. The third planet in the system has a 
period of 46.78 ± 0.16 days and a minimum mass of 8.7 ±1.1 Mg, at 0.234±0.002 AU from the star. Its eccentricity is 0.32 ± 0.14. 

The period of this planet is close to the rotational period of the star estimated from variations of activity indicators (42.3 ±0.1 days). 

The planetary origin of the signal is, however, the preferred solution as no indication of variation at the corresponding frequency is 
observed for activity-sensitive parameters. Finally, a fourth additional longer-period planet of mass of 62 ± 6 Mg orbits the star in 
1190 days, on an eccentric orbit (e = 0.27 ± 0.11) at a distance of 2.14±0.27 AU. 

Key words, planetary systems: super-Earth - techniques; radial velocity - techniques: photometry - stars:individual: HD 219134 - 
binaries: eclipsing - instrument: HARPS-N 


1. Introduction 

The statistical analysis and occurrence rate of the small-mass 
planets in the range of super-Earth to Neptune are discussed in 

* The photometric time series and radial velocities used in this work 
are available in electronic form at the CDS via anonymous ftp to 
cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi- 
bin/qcat?J/A±A/ 


several studies with the main motivation of better understand¬ 
ing the orbital and physical characteristics of this population of 
planets, in order to constrain their formation processes ( |Mor- 
Idasini et al.||2012[ |Benz et al.|2014| l. The radial-velocity planet 
search program with the HARPS spectrograph on the ESO 3.6- 


m telescope (Mayor et al. 2003| 

201 lb, a similar survey with 

the Keck telescope (|Howard et al. 

201011, and the NASA Kepler 

transit space mission (Borucki et a 

. 201 lb have in particular con- 
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tributed in a tremendous way to our knowledge of the population 
of small-mass/size planets around solar-type stars. 

The most common planets detected by the Kepler mission 
peak around 2 R@ ( [Howard et al.|2012)|Fressin et al.|2013HMarcy| 
let al.|2014| l. Even if such planets do not exist in our Solar Sys¬ 
tem, they are found around more than 30 % of solar-type hosts. 
Furthermore, many of them are found in very coplanar multiple 
systems ( [Fabrycky et al.|20T4{ [Howard et al.|2012[|Figueira et al?] 
2012^, tightly packed close to the central star, a new challenge 


to explain for planet formation models ( [Ogihara et al.[[20T5] l. 
Over the last 12 years the extra-solar planet zoo has also been 
supplied with low-mass planets detected by the HARPS GTO 
planet search program and successive subsequent ESO Farge 


programs. Notable examples include HD 40307 ([Mayor et al. 
[2009ll, HD 10180 ( [Lovis et al.|20TTbll, HD 20794 and HD 85512 


( Pepe etal.[2011[ l, or Alpha CenB ( Dumusque et al.[2012[ l. Anal¬ 
ysis of the results by Mayor et al.[ ( 201 l[ l provided a list of addi¬ 
tional super-Earth and mini-Neptunes unveiled by the survey, as 
well as a first statistical analysis of the properties and occurrence 
rate of the super-Earths and Neptune-mass planets around solar- 
type stars. These preliminary findings were confirmed by the 
fantastic statistics and precision of the Kepler detections (|Fressin[ 
let al.|2013» . 

After a decade focusing mainly on the detection and the de¬ 
termination of occurrence rate of low-mass planets, a significant 
observational effort is now dedicated to the planet characteriza¬ 
tion. Transit results provide the planet radius and in combination 
with radial velocities, the mean density of the planet can be de¬ 
rived. Spectral features of exoplanet atmospheres may also be 
revealed by space and ground-based high-resolution transmis¬ 
sion spectroscopy in the visible and near infrared. One of the 
main difficulties of such characterization is the availability of 
bright targets favourable for follow-up observations. In partic¬ 
ular exquisite planetary physical parameter determination will 
be key in lifting the intrinsic degeneracy in the determination 
of the planet composition when several components (gas, sili¬ 
cates, metals) and chemical species are mixed in the planet in¬ 
terior. This was a strong driver for transit search follow-up of 
super-Earths detected with radial velocities around very bright 
stars, from space with the MOST, Hubble or Spitzer space tele¬ 
scopes ( Gillon et al.|2012 2010[ Demory et al.|2015a) l, and for 
the development of a new generation of space transit missions to 
be launched within the coming decade: CHEOPS/ESA ([Fortier[ 
let al.|2014|l, TESS/ NASA ( [Ricker et al.[M4| , and PLATO/ESA 
( Rauer et al.|2014[ l. 

In order to start to fulfill the need for good estimates of planet 
physical parameters and considering the successful achieve¬ 
ments with HARPS, a similar design has been implemented for 
the HARPS-N spectrograph. The instrument was built to become 
an efficient exoplanet hunter in the northern hemisphere, specif¬ 
ically aimed at the follow-up and characterization of the bright¬ 
est transit candidates from the Kepler mission. HARPS-N is a 
high-precision echelle spectrograph located at the 3.6 m Tele- 
scopio Nazionale Galileo (TNG) at the Roque de los Muchachos 
Observatory on La Palma, Spain. It was built by an international 
consortium of Swiss, Italian, UK and US partners, led by Geneva 
university. It began operations in August 2012. Eighty nights per 
year for five years were granted to the consortium for a science 
program including the follow-up of Kepler candidates at high 
precision, and an additional survey, the Rocky Planet Search 
(RPS), to search for small planets around bright and quiet nearby 
stars in the northern hemisphere when Kepler field is down. 

While waiting for the new generation of ground- and space- 
based transit searches targeting bright stars, planet density esti- 
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mates from the follow-up of the Kepler candidates with HARPS- 


N have already been presented in several papers ([Pepe et al. 

2013 

Dumusque et al.[2014[[Bonomo et al.[2014[ Dressing et al. 

2015 

1 . They are populating the low-radius regime of the mass- 


radius relation for small exoplanets with candidates for which 
precise mass and radius could be determined. 


We present in this paper the first results from the Rocky 
Planet Search: a planetary system around HD 219134, composed 
of 3 inner super-Earths and an outer sub-Saturn mass planet, with 
the 4 planets spread between 0.038 and 2 AU from the central 
star. Thanks to the high precision photometric observations from 
the Spitzer space telescope, the inner planet is furthermore ob¬ 
served to transit in front of the star. After a short presentation 
of the Rocky Planet Search program in Sect.j^ and the stellar 
properties of HD 219134 in Sect.j^ the spectroscopic and pho¬ 
tometric observations of the new system are described in Sect.|^ 
an d0 Sect. [^provides a discussion of the stability of the system 
and of the composition of the inner planet. We finally draw some 
conclusions in Sect.jT] 

2. The HARPS-N Rocky Planet Search program 

The GTO (Guaranteed Time Observation) granted to the 
HARPS-N consortium is dedicated to two programs: i) the con¬ 
firmation of the planetary nature and the characterization of Ke¬ 
pler candidates and ii) a Rocky Planet Search (RPS). The aim of 
the RPS program is to perform a systematic search for low-mass 
planets around nearby quiet stars visible from the northern hemi¬ 
sphere, through an intensive monitoring of the radial velocity of 
quiet stars at very high precision. HARPS-N is stabilized and 
well-controlled in pressure and temperature to minimize instru¬ 
mental drifts, and so to secure sub-m/s radial velocities over long 
periods of time. More details on technical aspects are presented 
in [Cosentino et al!] ( [2012[[2014| l. 

The sample: The first step in establishing the RPS program 
was to select a sample of stars best suited for long-term pre¬ 
cise radial-velocity measurements. Uncertainties in such obser¬ 
vations are mainly linked to noise from three different origins: 
photon noise, instrumental noise, and stellar intrinsic signals. 
The first important criterion for this program was thus to concen¬ 
trate on bright stars, in our closest neighbourhood. This is also 
in the interest of potential follow-up studies for characterization 
of the planet properties. The sample was thus selected from the 
Hipparcos catalogue within a distance limit set to 18 P^] The 
second step was to focus on quiet stars. Based on the CORAVEL 
data and investigation in the literature, we rejected known spec¬ 
troscopic binaries as well as stars with V sin / > 4.5kms“^ and 
stars with log > -4.69 from our sample. In addition, K 

dwarfs are favoured in our sample compared to G dwarfs be¬ 
cause of their lower level of stellar "noise" (p-mode, granulation 
and activity jitter; [Dumusque et al.[201 Ibj i. Their habitable zone 
is also closer in. Adding all the criteria together, we ended up 
with a sample of 51 quiet stars with a range of spectral types 
from G8 to MO. They are listed in Table[T| 

Observational strategy: To minimize the effect of stellar 
noise with short typical time scales (p-modes and granula¬ 
tions), we applied the observational strategy implemented for 
the HARPS very high-precision observations ( [Pepe et al.|2011[ l. 
With this strategy, we observe each target with 15 minutes expo¬ 
sure time to damp the p-modes effect and make a 2nd (and possi- 

' The limit was chosen iteratively in order to have a sample large 
enough and covering a full range of right ascension, with an over den¬ 
sity during winter when the Kepler field in not visible 
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Target ID 

Distance [pc] 

V 

Spec.Type 

HD38 

11.75+0.38 

8.20 

K2 

HIP1368 

14.99+0.25 

8.99 

MO 

HD3651 

11.11+0.09 

5.88 

KOV 

HD4628 

7.46+0.05 

5.74 

K2V 

HD 10476 

7.47+0.05 

5.24 

KIV 

HD 10436 

13.43+0.20 

8.42 

K5V 

HD16160 

7.21+0.05 

5.79 

K3V 

HD 19305 

14.77+0.39 

9.07 

MOV 

HD232979 

10.19+0.12 

8.62 

K8V 

HD32147 

8.81+0.06 

6.22 

K3V 

HIP27188 

13.69+0.21 

9.02 

MO 

HD41593 

15.45+0.16 

6.76 

KO 

HD47752 

17.52+0.35 

8.08 

K2 

HD48948 

16.40+0.38 

8.59 

MO 

HIP36357 

17.55+0.29 

7.73 

K2V 

HD62613 

17.04+0.11 

6.55 

G8V 

HD65277 

17.46+0.39 

8.05 

K5V 

HD65583 

16.80+0.16 

6.97 

G8V 

HIP39826 

17.12+0.54 

9.41 

MO 

HIP42220 

13.90+0.24 

9.28 

M2 

HIP43534 

16.49+0.37 

9.26 

K5 

HD79210 

6.19+0.20 

7.64 

MOV 

HD79211 

6.27+0.26 

7.70 

K2 

HD79969 

17.53+0.23 

7.20 

K3V 

HD84035 

17.78+0.31 

8.13 

K5V 

HD88230 

4.87+0.20 

6.60 

K8V 

HIP51525 

15.67+0.27 

8.85 

K7V 

HIP51547 

17.49+0.42 

9.63 

MO 

HD94765 

17.55+0.26 

7.37 

KO 

HD97101 

11.93+0.15 

8.31 

K8V 

HD97503 

17.95+0.41 

8.70 

K5V 

HD99492 

17.99+0.47 

7.58 

K2V 

HD 103095 

9.16+0.07 

6.42 

G8V 

HD110315 

14.30+0.22 

7.91 

K2 

HDl 11631 

10.78+0.11 

8.49 

M0.5V 

HD 122064 

10.10+0.06 

6.49 

K3V 

HD128165 

13.42+0.12 

7.24 

K3V 

HD 144579 

14.37+0.12 

6.66 

G8V 

HD 147379 

10.66+0.14 

8.61 

MOV 

HD151288 

9.77+0.08 

8.10 

K7V 

HD158633 

12.80+0.08 

6.44 

KOV 

HD 157881 

7.72+0.06 

7.54 

K7V 

HD 166620 

11.10+0.07 

6.38 

K2V 

HD173818 

14.12+0.25 

8.81 

K5 

HD 185144 

5.77+0.02 

4.67 

KOV 

HD 184489 

14.47+0.35 

9.35 

K5 

HD 190007 

13.11+0.18 

7.46 

K4V 

HD200779 

14.84+0.30 

8.27 

K5 

HD201091 

3.48+0.02 

5.20 

K5V 

HD201092 

3.50+0.01 

6.05 

K7V 

HD219134 

6.53+0.03 

5.57 

K3V 


Table 1; List of the RPS targets in the GTO program on HARPS- 
N drawn from the Hipparcos catalogue and van Leeuwen (2007|l. 


bly a 3rd) measurement of the target well-spread in time during 
the night to damp the effect of granulation (for more details see 
[Dumusque et al.|2011b| l. For bright stars, in order to avoid sat¬ 
uration the 15 minutes on target are split into sub-observations. 
Typical average signal-to-noise ratio (SNR) of the spectra ob¬ 
tained are between 200 and 400, measured at T = 550 nm. 



RV dispersion [m/s] 


Fig. 1; Histogram of the dispersion of the radial-velocity mea¬ 
surements of stars in the RPS program on HARPS-N. 


With continuous monitoring during the past two years of op¬ 
eration of HARPS-N we have acquired hundreds of precise ra¬ 
dial velocities of the stars in the RPS sample. The typical photon- 
noise precision per [sub-]observation is about 0.3 ms *. In Fig.[^ 
we show the radial-velocity rms for all the stars in the RPS pro¬ 
gram. This rms naturally includes photon noise, instrumental ef¬ 
fects (telescope, spectrograph, detector), stellar intrinsic "noise", 
and of course signatures of still undetected planets. With a mode 
around 1.48 ms ', the distribution is very similar to the one ob¬ 
tained for the sample of solar-type stars followed at high preci¬ 
sion with HARPS in the southern hemisphere. 

We present here the first result of the HARPS-N RPS pro¬ 
gram; the discovery of a planetary system with a transiting planet 
around HD 219134. 


3. Stellar characteristics of HD 219134 


The very bright star HD 219134 has been extensively studied 
in the literature. Basic photometric properties of this star were 
taken from the Hipparcos catalogue ( |van Leeuwen|2007] ). A pre¬ 
cise estimate of the radius of the star, R = 0.778 + 0.005 Rq, 


is available from interferometric measurements (Boyajian et al. 
|2012| l. Combined with a luminosity of 0.265 + 0.002 Lq, the 
Stefan-Boltzmann law gives T^ff - 4699 + 16K. We will adopt 
this value for our photometric analysis of the transiting planet in 
the system in Sect.|^ 

Over the past few years, the star has been the object of sev¬ 
eral studies aiming at determining photospheric parameters and 


chemical abundance analyses (|Valenti & Fischer|2005| Mishen-| 

ina et al.|2008 

2012[|Kovtyukh et al.|2003| Ramirez et al.|2013 

Prugniel et al. 

201111. Results of these studies are reported in Ta- 
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Technique 

Teff [K] 

log(g) 


[Fe/H] 

SME* 

4835 + 44 

4.56 + 0.06 


0.09 + 0.03 

LDP2 

4889 + 20 

4.60 + 0.20 


0.10 + 0.05 

Photometry^ 

4833 + 60 

4.59 + 0.02 

0.26 

0.00 + 0.06 

ULySS^ 

4715+44 

4.57 + 0.07 


0.06 + 0.04 

EWs^ 

4820 + 61 

4.62 + 0.17 

0.35 

0.12 + 0.04 

SPC^ 

4941 + 50 

4.63 + 0.10 


0.11+0.08 


Valenti &■ Fischer {2005\; ‘ Mishenina et al.\{2008\; Ramirez et al.\ 


( 2013) ; ‘ Prugniel et al. [2011 1; ^Thispaper 


Table 2; Comparative summary table for the atmospheric param¬ 
eters derived for HD 219134. 


We also derived the atmospheric stellar parameters directly 
from HARPS-N spectra using two approaches; one based on 
Equivalent Width (EW) determination and one using the Stellar 
Parameter Classification (SPC) tool. 

Eor the Equivalent Width (EW) approach, we followed the 
procedure described in Sect. 3 of |Dumusque et aT ( 2014| l. We 
used the 2014 version of the local thermodynamic equilibrium 
code MOOG ( [Sneden |1973|l together wi th the Kurucz grid of 
atmosphere models (Castelli & Kurucz 2004[ [Kurucz 1992|l, 
while oscillator strength values in the line list from [Sousa et al.j 
( [201 l[ l were updated in accordance with the solar iron abun¬ 
dance loge(Eel) = 7.50 from Asplund et al.[(2009i. In order 


to construct high-qualitity spectra, spectra within the same visit 
(sub-observations) were co-added, obtaining at the time of the 
analysis 88 spectra with SNRs between 400 and 700. EWs were 
then measured for each of them and their mean and rms were 
used for the atmospheric parameter determination. We only re¬ 
tained the lines within the range 5 mA < EW < 100 mA and 
with a dispersion lower than either 1 mA or 5 % of the mean 
EW. Despite the quality of our spectra, the preliminary analysis 
resulted in a poor determination of the microtubulent velocity 


Eollowing the calibration from Ramirez et al. ( [2013 1 and Tsan- 
[taki et alT| ( [2013| l, we decided to fix the microturbulent velocity 
to = 0.35. The resulting atmospheric parameters are reported 
in Table with the gravity from Eell lines already increased by 
A log g - 0.39 + 0.15 according to the calibration in Sect. 3.1 of 
[Mortier et al.| ( [201^ . Note that the derived parameters are domi¬ 
nated by systematic errors (e.g. choice of the oscillator strength, 
continuum placement) rather than random errors on EW mea¬ 
surements. 

Einally, we used the SPC tool, an alternative approach to de¬ 
rive atmospheric parameters by cross-correlating an observed 
spectrum with a library of synthetic spectra (Buchhave et al. 
|2012[[20T4l l. With SPC we obtain an effective temperature Teft = 
4941 + 50 K, a surface gravity of 4.63 + 0.10 and V sin I = 
0.4 + 0.5 kms '. The metallicity derived by SPC is 0.11 + 0.08, 
from a mix of metallic absorption lines in the wavelength range 
between 5050 to 5360 A. The derived values are reported in Ta- 
ble|2|as well. 

The mass of HD 219134, M* = 0.78 + 0.02 Mq, was es¬ 
timated through the Synthetic Clusters Isochrones & Stellar 
Tracks tool (SYCLISTf] SYCLIST allows the user to deter¬ 
mine stellar parameters through a Bayesian-based interpolation 
of the grid of Geneva stellar evolution models. Eor HD 219134, 
we used the Tetj taken from Boyajian et al.[ ( 2012| l, the appar¬ 
ent magnitude and parallax from the Hipparcos catalogue dvanl 


Leeuwen [2007 [l, and the metallicity derived from our spectral 


http://obswww.unige.ch/Recherche/evoldb/index/ 
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Parameter 



Ref. 

V 


5.57 

1 

B-V 


0.99 

1 

SpTyp 


K3V 

1 

n 

[mas] 

152.76 + 0.29 

1 

L* 

[Lo] 

0.265 + 0.002 

2 

R* 

[Ro] 

0.778 + 0.005 

2 

Teff 

[K] 

4699 + 16 

2 

log g 

[dex] 

4.63+0.10 

this paper 

V sin i 

[kms-'] 

0.4 + 0.5 

this paper ** 


[kms-'] 

0.35+0.19 

this paper * 

[Fe/H] 

[dex] 

0.11 +0.04 

this paper 

M* 

[Mo] 

0.78 + 0.02 

this paper, 2 

<^os(R'hk> 


-5.02 + 0.06 

this paper 

Prot 

[days] 

42.3 + 0.1 

this paper 


Table 3: Stellar paramete rs taken from * the Hipparcos cata- 
logue (van Leeuwen|2007 1 , ^interferometric measurement (Boy 


[ajian et al.|2012| l and estimated from HARPS-N spectra adopting 
*EWs and “SPC and stellar evolution models SYCLIST for M* 
(see text). 


analysis as input parameters. Eor the latter we used the average 
of the values obtained from the EWs and SPC approaches. 

The mean activity index log R^^ is estimated from the 
HARPS-N spectra. We derive an average value of -5.02 with 
a dispersion of 0.06 dex. To estimate the rotational period of the 
star we used a periodogram analysis of the activity indicators 
(log CCE EWHM and CCE bisector span time series) that 
yields a most significant peak at 42.3 days (see Sect. |4.2.^ . 

Table [^summarises the stellar parameters obtained from the 
various catalogues and analyses mentioned above, with the final 
values selected for the star when several estimates are available. 


4. Spectro-velocimetric observations 

4.1. The HARPS-N data 

HD 219134 has been monitored by the HARPS-N spectrograph 
for close to three years, from BJD= 2456148.7 (9th of Aug. 
2012) to BID = 2457195.7 (22nd of June 2015). To mitigate the 
effects of stellar oscillations, the strategy was to stay on target for 
10 to 15 minutes. Because of the brightness of the star, the obser¬ 
vations were split in several individual sub-observations to avoid 
saturation. We obtained a total of 481 data points spread over 99 
epochs (nightly averaged values). The nightly averaged values 
are displayed in Eig.j^ Observations were performed using the 
simultaneous thorium calibration technique. The mean uncer¬ 
tainty on the individual RVs due to photon noise and known cal¬ 
ibration noise is 0.4 ms '. This corresponds to an average SNR 
of 268 per pixel at A- 550 nm. The data reduction was carried 
out with the latest version of the HARPS-N pipeline (DRS 3.7) 
using the cross-correlation technique ([Pepe et al.|2002||Baranne| 
et al.|1996| l. On top of precise radial velocities, the pipeline di¬ 
rectly provides parameters estimated from the cross-correlation 
function (CCE) of the spectrum: CCE full width at half maxi¬ 
mum (EWHM), CCE contrast (minimum vs continuum), CCE 
bisector span inverse slope, and Call activity index S and log 

Rntr (Fig-El- 

The raw rms dispersion of the radial velocities is 3.57 ms *. 
Once de-trended from an obvious drift, the observed radial ve¬ 
locities still show a dispersion of 2.77ms *, which is signifi- 
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JD - 2450000.0 [days] 

Fig. 2: Radial-velocity (top) and activity index (bottom) time se¬ 
ries for HD 219134. 


candy above the typical dispersion of such quiet stars, calling 
for a search for additional coherent signals in the data. 

4.2. Data analysis 
4.2.1. Data Modelling 

The first step of the radial-velocity data analysis consists in 
identifying significant periodic signals in the data. This was 
done using the General Lomb-Scargle periodogram algorithm 
(GLS, |Zechmeister & Kurster|[2009] l applied to the nightly av¬ 
eraged radial-velocity measurements to which a systematic er¬ 
ror of 1 ms ' was quadratically added. False alarm probabilities 
were estimated through a bootstrap approach by permuting the 
nightly averaged data. Once a significant peak was located at 
a given period, the corresponding Keplerian was adjusted and 
removed. The process was repeated several times until no sig¬ 
nificant peak remained. For multiple Keplerians, all parameters 
were re-adjusted at each step of the analysis. During the analy¬ 
sis, we identified one radial velocity outlier (5 cr at jdb=56829.7) 
and decided to remove it to allow for a robust frequency analysis. 

Following this procedure, a first low frequency oscillation at 
1190 day was identified with a semi-amplitude of 4.5 ms“' and 
a 0.8 % false alarm probability (FAP). After removing the corre¬ 
sponding Keplerians, two highly significant peaks were seen in 
the periodogram, at 3.09 and 46.78 days, with semi-amplitudes 
of 2.33 and 1.94ms ', and FAPs smaller than 0.01 % and 1 %, 
respectively (Fig.[^. One additional peak remains at 6.76 d with 
a significant FAP smaller than 1 % and a semi-amplitude of 
1.1 ms '. No more periodic signals remain in the data after sub¬ 
traction of the corresponding 4 Keplerian model. This purely fre- 
quentist approach, using the nightly average data, has the advan¬ 
tage of being "simple and fast" and is also conservative in terms 
of detection limits. 



5 . 10.0 20 . 50 . 100.0 200 . 500 . 1000 . 

Period [days] 


Fig. 3: Generalized Lomb Scargle periodogram of the radial ve¬ 
locities of HD 219134, after removing - from top to bottom - 
a long period Keplerian, and successively the 3.09 d signal, the 
46.78 d signal and finally the 6.75 day Keplerian. Dashed lines 
indicate 0.1%, 1%, and 10% level of false alarm probabilities. 


To double check the purely frequentist approach, we also 
analysed models of different complexity ranging from 1 to 4 Ke¬ 
plerian. The best model was chosen by comparing the Bayesian 
Information Criterion (BIC, Kass & Raftery|1995| l between the 
different models; 


BIC 2 log T-H Aparam Ittg Nmeas » 

where logL is the maximum of the log likelihood, Aparam the 
number of free parameters in the model, and A^eas is the number 
of measurements. In model comparisons based on the Bayes fac¬ 
tor, when the prior information is smaller than the information 
provided by the data, the BIC (a.k.a Schwarz criterion) indicates 
that the model with the highest probability is the one with the 
minimum BIC. A summary of the comparison between a few 
models is provided in Table]^ The 4 Keplerian model has the 
lowest BIC and differs from the 3 Keplerian model by ABIC=34. 
It is usually assumed that a ABIC of 20 between two models is 
considered as strong evidence in favour of the most complex 
one which leads us to adopt the 4 Keplerians as our best solution. 

To obtain robust confidence intervals for the Keplerian pa¬ 
rameters as well as an estimate of the additional noise present 
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Model (Period [days]) 

BIC 

ABIC 

xi 

^o-c 

K1 (1190) 

611 


6.35 

2.63 

K2(1190, 3.09) 

389 

221 

3.90 

2.00 

K3 (1190, 3.09, 46.8) 

275 

114 

2.47 

1.54 

K4(1190, 3.09, 46.8,6.76) 

241 

34 

1.89 

1.31 


Table 4; Models tested on the 98 nightly averaged radial veloci¬ 
ties with an additional noise of 1 ms ' quadratically added to the 
radial-velocity error bars. The 4 Keplerian model is preferred 
since it has the lowest BIC, ^nd residual dispersion. These 
values are estimated based on the maximum likelihood. 


in the data (nuisance parameter below), we further probe the 
parameter space with a Markov Chain Monte Carlo algo¬ 
rithm (MCMC) with Metropolis-Hasting. An abundant litera¬ 
ture discusses in much detail the implementation of MCMC 
posterior sampling (Andrieu & Thoms ( 2008)1 for a pure sta- 
tistical approach; |Gregory| (|2005a|b|), |Collier Cameron et ah 


( |2007| l, and |Pollacco et ah ~ ^2008| l for exoplanet searches). 
Our MCMC prohes the following set of parameters; log p, 
Vecoscu, V^sinm, log/f, To (the mean longitude at a given 
epoch) while the noise model follows a simple normal law with 
standard deviation derived from the observation errors and a nui¬ 
sance parameter + ■*^)- Jeffrey’s priors are used for the 

period, the radial-velocity semi-amplitude, and the nuisance pa¬ 
rameter while uniform priors are used for the other parameters. 


4.2.2. Orbital solution 

The orbital elements corresponding to the four Keplerian model 
and transit timing constraint are listed in table while the phase 
folded radial velocities are displayed in figures and The 
shortest period signal, with a radial velocity semi-amplitude of 
2.33 ms“' at 3.09 days, is clearly identified as a super-Earth 
(4.46 M®) in a circular orbit. One additional low amplitude sig¬ 
nal (K=l.lms ') is present at 6.76 days, corresponding to a 
planet with minimum mass of 2.67M®. The existence of this 
signal is corroborated by both the FAP and the BIC estimators 
but its strength might be affected by the un-perfect modelling of 
the long-period signal (see below). 

The third signal with a period of 46.78 days and a semi¬ 
amplitude of 1.94 ms ' corresponds to a super-Earth of 8.67 M® 
with an eccentric orbit of e = 0.32. We show, in the following 
section, that this period, which is close to the rotational period 
of the star (42.3 days), is not of stellar origin. Its relatively high 
eccentricity may also originate from an un-perfect modelling of 
the outer signal. 

The long period radial velocity trend clearly shows two ex¬ 
trema, which favour a Keplerian instead of a polynomial drift, 
with a period of 1190 days and a semi-amplitude of 4.5 ms“^. 
As explained below, its origin cannot be linked to any long term 
stellar activity fluctuation and is undoubtedly of planetary origin. 
Due to the incomplete coverage of the orbit, both the eccentricity 
and the period of the planet remain moderately constrained but 
still leads to a good estimate of its minimum-mass, i.e. 62 M®. 


4.2.3. Periodic signals: planet vs activity-related origin 

Long-period magnetic cycle: The radial-velocity and activity 
index log 7?)^ „ measurements of the star HD 219134 are dis¬ 
played in Eig.]^ Although both exhibit long-term variations, the 
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Eig. 4; Phase-folded radial-velocity measurements of 
HD 219134 with the corresponding Keplerian model (solid 
line) for each of the 3 inner super-Earths, after removing the 
contribution of all the other planets in the system. Erom top to 
bottom, we have the 3.09-, 6.76-, and 46.78-day periods. 


log time series is very stable over the hrst 750 days of the 
survey and slightly increases over the last 300 days whereas the 
radial velocity constantly decrease over the 1050 day observa¬ 
tion time span. No correlation is observed between RVs and log 
as shown on Eig.l^ The long term trend of the log 
is most probably relatea to the magnetic cycle of the star (|L,o- 


vis et al. 1)201 la) Dumusque et al.))2011a Meunier & Lagrange 

2013) 1 while the observed 1190 day period corresponds to a long 
period companion. Moreover, at the observed level of activity, 
we expect a low impact of stellar activity on the observed radial 
velocities. 
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Model 


HD219134b 

K4+ Af(0, 
HD219134C 

f + ^2) 

HD219134d 

HD219134e 

P 

[days] 

3.0937 + 0.0004 

6.765 + 0.005 

46.78 + 0.16 


K 

[m/s] 

2.33 + 0.24 

1.09 + 0.26 

1.94 + 0.29 

4.46 + 0.52 

do 

[deg] 

82 + 8 

295 + 20 

98+16 

206/2^ 

T, 

[BJD-2400000] 

57126.7001 +0.001 

57129.46 + 0.45 



\[e. cos(w) 


0.05 + 0.19 

0.17 + 0.26 

-0.43 + 0.18 

0.21 + 0.23 

Ve. sin(w) 


-0.11+0.21 

-0.03 + 0.31 

0.03 + 0.21 

-0.35 + 0.24 

e 


Q"QQ+(rn 

^•^^-0.00 

q7v)T(06 

^•^^-0.00 

0.32 + 0.14 

0.27 + 0.11 

CO 


undehned 

undehned 

143 + 33 

288 + 45 

nipi sin i 

[Me] 

4.46 + 0.47 

2.67 + 0.59 

8.67+ 1.14 

62 + 6 

a 

[AU] 

0.0382 + 0.0003 

0.064 + 0.001 

0.234 + 0.002 

p 14+0.43 

y 



-18.4203 kms ^ +0.6 ms * 


^meas 



98 



S 

[m/s] 


1.18 + 0.06 



Table 5; Orbital solution for the 4 Keplerian model (K4) and planet inferred parameters for the system around HD 219134. T, is the 
expected date of transit, do is the mean longitude at the time 2457126.7001 day, corresponding to the transit timing (see Sect.|^ 



BJD - 2450000.0 [days] 


Fig. 5; Radial-velocity measurements as a function of time of 
HD 219134 with the corresponding Keplerian model (solid line) 
for the outer planet in the system, after removing the contribution 
of the 3 inner ones. 


Rotation period: Although different, the period of the third 
planet (46.8 days) is not very far from the rotational period of 
the star estimated from coherent variations of the activity indi¬ 
cators. A periodic variation around 42.3 days is indeed observed 
in the log as well as in the CCF FWHM and in the bisector 
span of the radial velocities, whereas no signal appears around 
46.8 days, as illustrated in Fig. [^showing the GLS periodograms 
of these parameters zoomed around the periods of interest. We 
can therefore safely consider the 42.3 days as a valid estimation 
of the rotational period of the star (Pro?). Conversely, no signal 
at 42.3 days appears in the radial-velocity data, supporting the 
planetary solution. 

In order to avoid any misinterpretation of the stellar activity 
as a planetary signal, we pushed our investigations a bit further. 
We hrst examined the yearly aliases ( |Dawson & Fabrycky|2010| l 
of the 46.8 and 42.3 days signals (Fig.|7jl and conhrmed that one 
period is not an alias of the other. In a second step, we consid¬ 
ered several subsets of the data to check the persistence of the 
planetary signal over time and to mitigate the possible effects of 
discontinuities in the data sampling. The same features are ob¬ 
served. This conhrms that the planetary signal at 46.8 days is 




6500 7000 


BJD-2450000 [day] 

Fig. 6; Plot of the star activity index as a function of the radial- 
velocity residuals after removing the 3 shorter-period planet con¬ 
tributions for HD 219134. No clear positive correlation is ob¬ 
served, contrary to what is expected in the case of magnetic cy¬ 
cles (Lovis et al. 201 la]l. 


present in the radial velocities at any time, and that in the same 
data no signal appears at the estimated (42.3 days). 

A hnal argument in favour of the planetary interpretation of 
the 46.8 day signal is provided in Fig.by the absence of a 
correlation between the logP^^ activity index and the residu¬ 
als around a 3-planet solution (leaving the 46.8 day period out). 
Such a correlation would be expected if the radial-velocity vari¬ 
ation is induced by activity-related spots or plages on the star 
surface. 


Article number, page 7 of 


13 




























A&A proofs: manuscript no. HD219134 



Fig. 7: From Top to Bottom the panels show the GLS peri- 
odograms of radial velocities, log CCF bisector span, and 
CCF FWHM of the 2-year data set, zoomed around the periods 
of interest. The thick red line in each panel, indicates the plan¬ 
etary period (46.8 days) and the thick green line indicates the 
stellar rotation period (42.3 days). The thinner lines are the cor¬ 
responding yearly aliases. 


6500 7000 

BJD-2450000 [day] 

Fig. 8; Radial velocities after subtracting the contributions of the 
planets at 3.09, 6.76 and 1190 days, displayed as a function of 
the logR^^ activity index detrended with a polynomial of de¬ 
gree 3. No correlation is observed supporting so the planetary 
interpretation of the 46.8-day signal still present in the data. 


5. Transit detection with Spitzer space teiescope 

5.1. Spitzer observations 

Thanks to its short-period of 3.09 d, HD 219134 b has an interest¬ 
ingly high geometric transit probability of ~9%. In Feb 2015, we 
requested 9.5hr of Director’s Discretionary Time (DDT) on the 
Spitzer Space Telescope to continuously monitor the 2-cr tran¬ 
sit window of the planet as derived from our analysis of the 
HARPS-N RVs. As demonstrated by its co-detection and sub¬ 
sequent studies of the transits and occultations of 55 Cnc e H 
|mory et al.|2011||2012[[2015b[|Gillon et al.|2012|l and its confir- 


mation of the transiting nature of HD 97658 b ( Van Grootel et al. 
2014|l, Spitzer is indeed a very powerful facility to search for and 


measure with high-precision the transits of low-mass RV plan¬ 
ets, thanks to its trailing orbit that allows monitoring the same 
star continuously for a complete transit window, and to its ultra- 
high photometric precision (a few dozens of ppm per time in¬ 
terval of 30 min for 55 Cnc). Our program was approved by the 
Spitzer Science Center (SSC) Director under the ID 11180 (PI; 
M. Gillon), and the SSC managed to schedule it for 2015 Apr 
14th, a few days before the end of the current visibility window 
of the star. 

We observed HD 219134 at 4.5 pm with the Spiteer/IRAC 
detector ( [Fazio et al.||2004] l in subarray mode (32x32 pixels 
windowing of the detector), the extremely fast Fowler sam¬ 
pling (~0.01s) maximizing the duty cycle and SNR. No dither¬ 
ing pattern was applied to the telescope (continuous staring). 
HD 219434 being an extremely bright star for Spitzer, we used 
the shortest available integration time of 0.01s, allowing the 
counts to remain in the linear regime of the detector. The ob¬ 
servations were performed from 2015 Apr 14th 01h58to llhl8 
UT. We used the recently introduced PCRS peak-up mode ([Grill- j 
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[mair et al.|2012[ [Ingalls et al.|2014[ ) which was previously used 
by Ballard et al. ( [2014[ l to estimate the infrared transit depth of 
Kepler-93 b. This mode provides enhanced accuracy in the posi¬ 
tion of the target on the detector, leading to a significant decrease 
of the so-called ‘pixel phase effect’ that is the most important 
source of correlated noise in high-SNR staring mode observa¬ 
tion with IRAC InSb arrays (e.g. [Knutson et al.|20()8| ). The run 
consisted of a 9hr-long science Astronomical Observational Re¬ 
quests (AOR) preceded by a short (30 min) AOR to allow the 
spacecraft to stabilize. 


5.2. Data reduction 


After download to Earth and basic calibration with the Spitzer 
pipeline S19.1.0, the images were made available to us by SSC 
through the Spitzer Heritage Archive (SHA) web interfac^ un¬ 
der the form of Basic Calibrated Data (BCD). Each subarray 
mode BCD is composed of a cube of 64 subarray images of 
32x32 pixels (pixel scale = 1.2 arc second). We used the fol¬ 
lowing strategy to reduce these BCDs. We first converted fluxes 
from the Spitzer units of specific intensity (MJy/sr) to photon 
counts, then aperture photometry was performe d on each sub - 


array image with the IRAF/DAOPHOlj^ software (Stetson 


1987 I. 


We tested different aperture radii, and selected 2.3 pixels as the 
radius minimizing at best the white and red noises in the resid- 


^ http://sha.ipac.caltech.edu 

IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in As¬ 
tronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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uals of a short data fitting analysis. The centre and width of the 
Point-Spread Functions (PSF) were measured by fitting a 2D- 
Gaussian profile on each image. The x - y distribution of the 
measurements was then looked at, and measurements having a 
visually discrepant position relative to the bulk of the data were 
then discarded. For each block of 64 subarray images, we then 
discarded the discrepant values for the measurements of flux, 
background, x- and y-positions using a lO-cr median clipping 
for the four parameters, and the resulting values were averaged, 
the photometric errors being taken as the errors on the average 
flux measurements. Finally, a 50-cr slipping median clipping was 
used on the resulting light curves to discard outliers (due, e.g., to 
cosmic hits). 

Our resulting light curve counted 9396 measurements. It is 
shown in Fig. with the evolution of relevant external parame¬ 
ters (PSF X- and y-center and PSF width, background). Its time 
sampling ~ 3.4i being much shorter than the structures of the 
expected transit and of the typical timescale of the Spitzer sys- 
tematics, we binned the light curve to time intervals of 30s for the 
sake of computational speed of the data analysis. Nevertheless, 
we verified with a shorter version of the data analysis procedure 
described below that our results are insensitive to the binning of 
the photometry. 



7126.6 7126.7 7126.8 7126 . 

BJD - 2450000. 


7126.6 7126.7 7126.8 7126.9 7127 

BJD - 2450000. 


5 . 3 . Data analysis 

We analysed the Spifzer photometric time-series with our adapta- 
tive MCMC code (see Gillon et al. 2012 and references therein). 
The assumed photometric model consisted of the eclipse model 
of [Mandel & Agol| ( |2002| l to represent the possible transit of 
HD 219134 b, multiplied by a baseline model aiming to represent 
the other astrophysical and instrumental effects at the source of 
photometric variations. We assumed a quadratic limb-darkening 
law for the star. We based the selection of the baseline model on 
the minimization of the Bayesian Information Criterion (BIC, 


|Schwarz|1978 

1 . 

Following 

Gillon et al.|(|2014]|, the instrumental models in 


eluded three types of low-order polynomials. The first one had as 
variables the x- and y-positions of the centre of the PSF to repre¬ 
sent the ‘pixel phase’ effect of the IRAC InSb arrays (e.g. |KnuF] 
[son et al.|2008| ). The second one had as variables the PSF widths 
in the x- and/or the y-direction, its inclusion in the baseline 
model strongly increasing the quality of the fit for Warm Spitzer 
photometry (see also Lanotte et al. 2014| l. The third, optional, 
function was a polynomial of the logarithm of time + a slope to 
represent a sharp decrease of the detector response at the start 
of the AOR (‘ramp’ effect, Knutson et al.|2008 1 . To improve the 
quality of the modelling of the pixel phase effect, especially the 
fitting of its highest frequency components, we supplemented the 
X- and y-polynomial with the Bi-Linearly-Interpolated Sub-pixel 
Sensitivity (BLISS) mapping method ( [Stevenson et al.|[2012| l. 
The sampling of the positions space was selected so that at least 
five measurements fall within the same sub-pixel. We refer the 
reader to jGillon et al.| ( |20T4) l for more details. 

Assuming no transit of HD 219134 b, the light curve cor¬ 
rected for the instrumental model described above showed a 
transit-like structure of ~ 300 ppm depth and lasting ~50 min 
(see Fig.[T0|), this structure remaining if higher order terms are 
used in the polynomial functions. Our first reflex was to check 
that this transit-like structure did not correspond to any odd 
behaviour of the external parameters, which was not the case 
(Fig.|g. Before identifying the structure to the searched transit, 
we nevertheless performed a large set of short MCMC analyses 
assuming different baseline models, and assuming or not that the 


Fig. 9; Evolution of the following measured parameters in the 
Spitzer images for HD 219134: PSF centre x- and y-positions 
(top), PSF X- and y-width (middle), background counts (bottom 
left), and stellar fluxes (bottom right). The times of the start and 
end of the detected transit are shown as red vertical lines. 


structure was a transit of the planet. For each baseline model, we 
computed the BIC difference between the best-fit models with 
and without transit to estimate the Bayes factor in favour of the 
transit hypothesis, multiplying the likelihood term 
9/91, the prior transit/no-transit probability ratio. In these tests, 
we multiplied the photometric errors by a Correction Factor (CF, 
see below) of 2.72, the highest value that we ever encountered in 
our past experiences of high-precision photometry with Warm 
Spitzer. It corresponds to an AOR targeting 55 Cnc (Demory et 
al., in prep.) for which the PCRS peak-up mode did not work 
properly. Doing so, we thus assumed that the HD 219134 light 
curve was strongly affected by some correlated noise poorly re¬ 
produced by our instrumental model. At the end, the lowest value 
that we obtained for the Bayes factor under this extreme hypoth¬ 
esis was ~ 1000 in favour of the transit hypothesis. We could thus 
conclude with certainty the transiting nature of the planet. 

We then performed a longer MCMC analysis to probe the 
posterior probability distribution of the transit parameters. The 
orbit of HD 219134b was assumed to be circular in this MCMC 
analysis. The jump parameters of the MCMC, i.e. the parameters 
randomly perturbed at each step of the Markov Chains, were the 
following. 

- The stellar mass M„, radius /?», effective temperature Tefr, 
and metallicity [Fe/H]. For these four parameters, normal 
prior probability distribution functions based on the values 
given in Table were assumed. 

- The planet/star area ratio dF - (RpjRff. 

- The impact parameter b - a cos ///?» of the planet, where a 
is the orbital semi-major axis and i is the orbital inclination. 
A full-transit configuration corresponds to 5 < 1 - RpfR,. 

- The time of inferior conjunction Tq for the planet, corre¬ 
sponding to the mid-time of the transit. 
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BJD - 2450000 


Fig. 10; Spitzer photometry divided by two different best-fit 
baseline models. The first one {top) aims to remove position- and 
PSF effects, and consists in a 4th order polynomial of the PSF x- 
and y-centres and widths, supplemented with the BLISS map¬ 
ping method ( [Stevenson et al.|[2012| l to remove high-frequency 
position effects. The second baseline model (bottom) adds to the 
first one a slope, and a quadratic function of the logarithm of time 
to model the sharp decrease of the counts at the beginning of the 
run (negative ramp). Both light curves are binned per 0.005d = 
7.2min intervals. 


The limb-darkening of the star was modeled by a quadratic law 
( |Claret||2()00[ ). Values for the two quadratic limb-darkening co¬ 
efficients Ml and M 2 were drawn at all steps of the MCMC from 
normal distributions with expectations and standard deviations 
drawn from the tables of jClaret & Bloemen] ( |201 Ij ) for the Spitzer 
4.5 yum bandpass and for the stellar atmospheric parameters 
drawn at the same step. 


Five chains of 100,000 steps were performed for each anal¬ 
ysis, their convergence being checked using the statistical test 
of jGelman & Rubin] ( 1992| l. They followed a preliminary chain 
of 100,000 steps performed to estimate the need to rescale the 
photometric errors, at the end of which the standard deviation of 
the residuals was compared to the mean photometric errors, and 
the resulting ratios I3„ were stored. represents the under- or 
overestimation of the white noise of each measurement. On its 
side, the red noise present in the light curve (i.e. the inability of 
our model to represent perfectly the data) was taken into account 
as described in jGillon et al.| ( |2()10j l, i.e. a scaling factor fir was 
determined from the standard deviations of the binned and un¬ 
bind residuals for different binning intervals ranging from 5 to 
120 minutes, the largest values being kept as y6,.. At the end, the 
error bars were multiplied by the correction factor CF - firXfi^.. 
The derived values for fir and /?„ were 1.30 and 1.01, resulting 
in CF = 1.32. 


Tab le|^ presents the resulting values plus error bars for the 
transit and planet’s parameters, while Fig 11 shows the light- 
curve corrected for the systematics, the best-fit transit model, 
and the residuals. 
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Fig. 11; Spitzer photometry divided by the best-fit baseline 
model and binned per 0.005d = 7.2min, with the best-fit transit 
model over-imposed in red. Below are shown the y-shifted resid¬ 
uals of the fit binned per 7.2 min and 30 min intervals. Their stan¬ 
dard deviations are, respectively, 57 ppm and 25 ppm. Over the 
light curve is shown the prior probability distribution function 
derived for the transit timing of HD 219134 b from our analysis 
of the HARPS-NRVs. 


Transit parameters 

Value 

Depth dF 

359 + 38 ppm 

Impact parameter b 

0.920 +0.010 R. 

Timing To 

2457126.7001 ± 0.0010 BJDtdb 

Duration W 

57.4 + 2.4 min 

Physical parameters 

Planet radius Rp 

1.606 +0.086 R® 

Orbital inclination i 

85.058 + 0.080 deg 


Table 6; Transit and physical parameters deduced from the 
MCMC analysis of the Spitzer photometry. 


6. Discussion 

6 . 1. Dynamical stability 

A very important and necessary a posteriori consistency check 
of a planetary system characterization is needed to verify that 
the dynamical evolution of the system is viable on the long 
term, assuring the persistence of the system from the end stage 
of its formation (when the protoplanetary disk disappears) till 
its observation today. Pure n-body integrations of the 4-planet 
system, using both the GENGA Code ( [Grimm & Stadel||2014| l 
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Fig. 12: Planet mass vs distance to the host stars for planets in 
close neighborhood. Transiting planets with good mass and ra¬ 
dius (density) determinations are shown as * in this diagram. 
The planets detected by radial velocity are also shown using 
the minimum-mass as a proxy for the mass (data taken from 
WWW. exoplane ts. org). 


launched through the DACE interface]^ and a fourth-order Her- 
mite scheme ( |Makino|1991| ), with initial conditions derived from 
parameters in Tablel^ and assuming coplanarity and an inclina¬ 
tion of 5 degrees from the transit observation, shows that the 
system is stable for more than 10^ orbits of the outermost planet. 

A more complete analysis, taking into account general rel¬ 
ativity and tides as well as longer-term secular effects, is be¬ 
yond the scope of this paper. As expected, a preliminary check 
indicates, however, that general relativity and tides will act to 
limit the secular growth of the eccentricities of the inner planets, 
favouring the long-term stability of the system. 


6.2. Bulk composition of the planet HD 219134 b 


HD 219134 b is the nearest transiting super-Earth known to¬ 
day. The radius and mass of the planet have been characterized 
to 6% and 9% accuracy, allowing us to place tight constraints 
on the bulk density. It thus will become one of the best targets 
for internal structure and atmosphere characterization with the 
Spitzer and HST and the up-coming space follow-up missions 
e.g. TESS, JWST, CHEOPS and PLATO. Its position in a stel¬ 
lar distance vs planetary mass diagram is provided in Eig.fT^ 
in comparison with the other most favourable known cases for 
characterization within 40 pc from the sun and up to 35 M®: 
GJ436b ([Butler et al.|2004||Gillon et al.|2007|), 55 Cnc e ([Gillon 


et al.|2012 |de Mooij k al.|2014|l, GJ 1214b (Charbonneau et al. 


2009)l, HD 97658 b poward et al.|[20TT] |Dragomir et al.||2012 


Van Grootel et al.|2014| l, GJ 3470 b ( Bonfils et al.|2012[ Demory 

et al.|2013| l and HAT-P-11 b ( [Dittmann et al.|20091 l. ' 

Combining the spectroscopic and space-based photometric 
data, we estimate a density of 5.89 + 1.17 gcm“^ for the planet. 
This density is consistent with the value of 6.90 gcm“^ that 
would be predicted for a 1.606 R© planet obeying the Earth-like 
compositional model presented in |Dressing et al. (|2015|l. 

The compositional tracks employed in Dressing et al. (2015 i 
(reproduced and updated in Eig.[T3]l are based on interior struc- 


^ DACE is a platform of the National Centre for Competence in Re¬ 
search ‘Planets’, that can be accessed at http:\\dace.unige.ch 


ture models by |Zeng & Sasselov] ( |2013| l that represent small 
planets as fully differentiated iron cores surrounded by lower 
density magnesium silicate mantles. These models provide a 
useful framework for comparing relative planet compositions, 
but the absolute core mass fractions are underestimated slightly 
because the Zeng & Sasselov ( |2013[ ) models do not incorporate 
the presence of lighter elements in the core and the inclusion 
of water in the mantle. Accordingly, the Earth-like composi¬ 
tional track presented in [Dressing et al. ( 2015) 1 corresponds to 
a model composition of 83% MgSi 03 and 17% Ee whereas the 


actual core mass fraction of the Earth is closer to 30%. Zeng 
et al. (2015|) have recently updated their model framework to ac¬ 


count for the presence of lighter elements in the core and the 
presence of water in the mantle. Employing the new models, 
we find that the population of highly irradiated dense planets 
(HD219134b, CoRoT-7b, Kepler-lOb, Kepler-36b, Kepler-78b, 
and Kepler-93b) are best described by a two-component iron- 
magnesium silicate model with a core mass fraction of approxi¬ 
mately 22-23%. 


7. Conclusion 

We have presented in this paper the first result from the Rocky 
Planet Search (RPS) program conducted with HARPS-N, as a 
planetary system of three inner super-Earths and one outer sub- 
Saturn planet, hosted by the bright and nearby quiet K dwarf 
HD 219134. The planet separations between 0.039 and 0.23 AU 
called for a search of potential transits of the inner planet(s) with 
the Spitzer space telescope. The successful detection of the tran¬ 
sit of HD 219134 b makes the star the closest (6.5 pc) and bright¬ 
est (V=5.5) star known to date with a transiting planet (super- 
Earth). This system is thus becoming one of the most favourable 
ones for follow-up observations, for better constraining the sys¬ 
tem architecture or aiming at the characterization of the planet 
physical properties. Eirst, it provides an exquisite constraining 
point in the mass-radius diagram (Eig.[T3]l for the composition 
of the planet, found to be of terrestrial-equivalent composition 
with a core mass fraction of the order of 22 %. The quality of 
the measurements of the radius, mass and then mean density ac¬ 
tually foreshadows what can be expected from the future transit 
missions in preparation that will target bright stars (CHEOPS, 
TESS, PLATO). We also know from Kepler results that multi- 
transiting systems of small-size planets are numerous. It is then 
now highly suitable to search for traces of transits of the other 
planets in the systems. Einally, even if a potential atmosphere 
around the planet is expected a priori to be tiny, the brightness 
of the system makes it worth trying to detect features of this at¬ 
mosphere in the UV, visible and NIR, from space and from the 
ground, especially in preparation for future measurements with 
larger facilities (JWST, TMT). 
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Fig. 13: Mass-radius relation for planets with radii smaller than 2.7 Re and with masses determined to a precision better than 20 % 
(updated from |Dressing et aT|p015| l). The shaded gray region in the lower right indicates planets with iron content exceeding the 
maximum value predicted from models of collisional stripping ( [Marcus et al.||201^ . The solid lines are theoretical mass-radius 
curves ( jZeng & Sasselov]|2013| l for planets with compositions of 100% H 2 O (blue), 25% MgSi 03 - 75% H 2 O (purple), 50% 
MgSi 03 - 50 % H 2 O (green), 100 % MgSi 03 (black), 50 % MgSi 03 - 50 % Fe (red), and 100 % Fe (orange). In this diagram, the 
position of HD 219134 b is almost overlapping the point for CoRoT-7b. It belongs to a group of planets including Kepler-36 b, 
Kepler-93 b, and Kepler-10 b. 
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